Non-equilibrium electromagnetic fluctuations: Heat transfer and interactions 
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The Casimir force between arbitrary objects in equilibrium is related to scattering from individual 
bodies. We extend this approach to heat transfer and Casimir forces in non- equilibrium cases 
where each body, and the environment, is at a different temperature. The formalism tracks the 
radiation from each body and its scatterings by the other objects. We discuss the radiation from 
a cylinder, emphasizing its polarized nature, and obtain the heat transfer between a sphere and 
a plate, demonstrating the validity of proximity transfer approximation at close separations and 
arbitrary temperatures. 
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The electromagnetic field in the space around bodies 
is stochastic due to quantum and thermal fluctuations. 
The basic formalism of Fluctuational Electrodynamics 
(FE), was set out over 60 years ago by Rytov and 
has been applied extensively since to diverse problems in 
radiative heat transfer 0, y[ and Casimir forces Q . FE 
starts with casting the current fluctuations in each body 
in terms of its dielectric properties, and proceeds to com- 
pute the resulting electromagnetic field. The improved 
precision of measurements of force and heat transfer at 
sub-micron scale have provided renewed incentive to ex- 
amine FE for objects at different temperatures In 
particular, when the size or separation of the objects is 
comparable to, or smaller than, the thermal wavelength 
(around 8 micron at room temperature), heat radiation 
and transfer will differ from the predictions of the Stcfan- 
Boltzmann law: The considerably larger near-field heat 
transfer, due to tunneling of evanescent waves, has been 
verified experimentally [a, 0] ■ Theoretical computations 
of heat transfer were only recently extended from two 
parallel plates Q or dipoles [U to two spheres @. The 
radiation of single spheres and plates has been studied 
by many authors [id For a cylinder, the emissivity 
restricted to waves traveling perpendicular to its axis has 
been addressed [11]. There are also recent computations 
of the non-equilibrium Casimir force between objects at 
different temperatures, for parallel plates 11311 . modulated 
plates 14| , as well as a plate and an atom [lfj . The limi- 



tation of these results to simple shapes and arrangements 
points out the need for approaches capable of handling 
more complex situations. 

Here, we derive a formalism for computing heat trans- 
fer and Casimir forces for arbitrary objects (compact or 
not) maintained at different temperatures. Generaliz- 
ing previous work on Casimir forces in equilibrium, our 
approach enables systematic description of FE of a col- 
lection of objects in terms of their individual scattering 
properties. For the non-equilibrium Casimir force, we can 
investigate interactions between compact objects where, 



unlike previous studies [13Hl5l| , the effect of a third tem- 
perature (of the environment) has to be taken into ac- 
count. In terms of new applications, we derive the heat 
radiation of a cylinder which is of interest for heated wires 
or carbon nanotubes [l6[. We also study the heat trans- 
fer between a sphere and a plate, the only geometry for 
which near field heat transfer has been measured [g, [7| . 

Consider an arrangement of N objects labelled as 
a = 1 . . . N, in vacuum at constant temperatures {T a }, 
and embedded in an environment at temperature T env . 
In this non-equilibrium stationary state, each object is 
assumed to be at local equilibrium with current fluctua- 
tions obeying the fluctuation-dissipation theorem (FDT). 
In the following, we derive the autocorrelation function 
C of the electric field E at frequency ui at points r and 
r' outside the objects, from which the Poynting vec- 
tor for heat transfer and the Maxwell stress tensor for 
Casimir forces can then be extracted. In equilibrium, 
with T a = T env — T, C is related to the imaginary part 
of the dyadic Green's function Gij by [H 
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where ax (to) 

tional to the occupation number of all oscillators of fre- 
quency u>, c is the speed of light, and h is Planck's con- 
stant. Zero point fluctuations which contribute ao(w) = 
- rt ( 4?r ) play no role in our discussion. We shall hence- 

= Gij(u>; r,r0. 
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forth employ the operator notation 



Since ImG = 
J2 a Ime a I = 



-GlmG" 
4lm(G- 



and using the identity 17 1 
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where e a is the com- 



plex dielectric function of object a and Go is the Green's 
function of free space, we obtain 

C eq (T) = Co + Y J C s a c (T) - orH-^GImGo J G*, 

' UJ 
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G« c (T) = a r (w)GIm £a G*, (2) 
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where Co = ao(cj)^ImG is the zero point term. The 
finite temperature contribution is thus a sum of TV + 1 
terms: Each C^(T) contains an implicit integral over 
sources within a and is identified with the field sourced 
by this object [l[; the scattering of this radiation by all 
other objects is accounted for by multiplying Ime Q on 
both sides with the full Green's function. The last term 
in Eq. ©, C env {T) = -ar^^GImG^G* , is hence 
identified with the contribution sourced by the environ- 
ment. 

A key assumption of FE is that in a non-equilibrium 
situation, the thermal current fluctuations inside each ob- 
ject are described by the FDT at the corresponding local 
temperature, and are independent of the impinging radi- 
ation from the other objects. Having identified the differ- 
ent sources in Eq. ([2]) , we can change their temperatures 
to arrive at the desired non-equilibrium generalization 

C neq (T env , {T a }) = C + Y, C s a c {T a ) + C env (T env ) 

a 

= C e "{T env ) + [C s a c {T a ) - C s a c {T env )] . (3) 

a 

The second form is obtained by considering the differ- 
ence of C neq (T env ,{T a }) from C eq {T env ) due to the de- 
viations of the object temperatures T a from T env . This 
form is useful because the equilibrium correlation can be 
regarded as known, and the number of sources is reduced 
from TV + 1 to N. Applying the formalism, e.g., to de- 
rive Casimir forces, the first term on the r.h.s. of Eq. Q 
yields the equilibrium force at temperature T env . 

The next step is to compute the radiation field of ob- 
ject a when isolated, i.e., before this field is scattered by 
the other objects, and with T env = 0. This is given by 
C a (T a ) = ax a (w)G ct Ime a G* where G Q is the Green's 
function of object a in isolation, and thus involves an 
implicit integration over the interior of object a. To 
employ multiple scattering techniques [l8j| . it is consid- 
erably more convenient to express C a (T a ) in terms of 
the T-operator or scattering amplitude T Q of the object. 
In equilibrium, the electric field correlator for the iso- 

2 

lated object C^ q (T a ) = aT„(u)^ImG a , contains radia- 
tion sourced (i) by the environment and (ii) by the object 
itself. The latter can be obtained by subtracting the con- 
tribution from the environment, which can be regarded 
as an additional material with e env — > 1, occupying the 
space complimentary to a \vi\ . Towards this calculation, 
we introduce a Green's function G Q with e a inside object 
a and e env outside, 
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Note that all sources for C^™"(T a ) are outside object a, 
and none of the Green's functions appearing in Eq. (|4]) 
contain points inside the object, which can thus be writ- 
ten in terms of T a as G Q = Go — GqT q Go [l8| (G q is a 



simple modification of G Q as a finite e env — 1 only changes 
the external speed of light). For computing the energy 
radiated by object a, one does not have to find C^ q (T a ): 
As a consequence of detailed balance, it does not con- 
tribute to the Poynting vector. 

Finally, to compute C^ c (T a ) in Eq. ([3]), we need to ac- 
count for scattering of the radiation emerging from a, by 
all other objects collectively designated by (3. Denoting 
their total T-operator by T^, by use of the Lippmann- 
Schwinger equation , we arrive at the final form 

C s a c {T a ) = O aJ C a {T a )Ol M , with (5) 

"^^-^ i-WV 

Expanding the resolvent leads to an alternating applica- 
tion of Tg and T a , corresponding to a sequence of scat- 
terings between the objects. Equations ([5]) and (j3)) con- 
stitute our non-equilibrium formalism. 

The correlator C neq enables computing the Poynting 
vector and Maxwell stress tensor, respectively given by 



S(r) 



c 

47T 



(E(w,r) x B*(w,r)) , 



where the arguments uj and r are omitted in the lower 
line. The heat H absorbed per unit time by object a, 
and the force F$ acting on this object in direction i, are 
then obtained by integrations of S and Tij over a surface 
o~ a enclosing only this object, as 

H a = — Re d> S • n Q da, F ijQ —Rcd> Tijn a j da, (6) 

where n a is the outward normal to the surface a a . 

As a first application we compare heat radiations from 
a single object, a plate, sphere or cylinder; the only 
shapes amenable to analytic treatment. As formulae for 
heat radiation of a sphere or a plate are available in the 
literature 0, [l(J , we focus on the cylinder where the cor- 
responding result is only discussed implicitly [ljjj]- For 
an infinitely long cylinder, T is represented in cylindrical 
wave functions [18j . indexed by (n, fc|| , P) where n is the 
multipolc order, fc|| the wave vector component along the 
cylinder, and P = E or P = M the polarization. The 
matrix element describes the relative amplitude of 

the scattered wave of mode (n, fc|| , P') emerging from an 
incoming wave (n, k^,P). We then find for the radiated 
heat of the single cylinder per length L, 



duj c 4 
aT{ujj- 
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P=E.M n=-oo 



J &||(Re[T^J + |T^J 2 + |Ti:g 2 ), (7) 
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FIG. 1. Heat radiation of a cylinder and a sphere of SiC>2, 
as function of R, normalized by the Stefan Boltzmann result, 
at T — 300K. The horizontal line shows the radiation of a 
SiC>2 plate. At and the smallest skin depth S in the relevant 
frequency range are marked on the 7?-axis. For the cylinder, 
the contributions of the different polarizations are shown. 



where P = M, if P = E and vice versa. In Fig. [TJ 
we compare the heat radiation of a plate (semi-infinite 
body), a sphere and a cylinder (both of radius R), all 
evaluated with optical data of silicon-dioxide (Si02), 
as used in experiments |(|. The radiation is normal- 
ized to the Stefan-Boltzmann law H = aT A A 1 with 
a = Tr 2 kg/(60h 3 c 2 ) and surface area A of the object. 

For thin cylinders and small spheres H is proportional 
to the volume, while in the opposite limit, it is propor- 
tional to the surface area, reflecting a finite skin depth 
(absorption length) S(u>) — c/^m^/euj): Thermal fluctu- 
ations at frequency uj within the object emit radiation 
which may be rc-absorbed on its way out. If S <C R, only 
thermal fluctuations near the surface lead to emerging 
radiation, while for 5 3> R, the entire volume contributes 
to H . An interesting feature of Fig. [1] is the intermediate 
range, where the sphere and cylinder emit more strongly 
than a plate of equal area, related to Mie resonances for 
the sphere [ll|. For R — » oo, i.e., when the wavelengths 
involved (roughly peaked around the thermal wavelength 
At = hc/ksT « 7.6/im), as well as skin depths, arc 
much smaller than the smallest dimension of the object, 
the classical (plate) result is approached. The asymptotic 
value is in these units denoted as cmissivity e(T) < 1. In- 
terestingly, the radiation from a cylinder is polarized, its 
parallel and perpendicular polarizations obtained from 
the P = E and P = M terms in Eq. ([7]), respectively. 
The predominant radiation of a thin cylinder is parallel 
and changes to perpendicular for R At- Both polar- 
izations become equal asymptotically as R — > oo. Polar- 
ization effects have indeed been observed experimentally 
for wires 12, 20J and carbon-nano-tubes |2l|, for which 



Now we consider multiple objects. To compute C neq 
involving spheres or cylinders, we also need to convert 
among bases appropriate to the different objects. For 
example, for the experimentally most relevant configura- 
tion of a sphere and a plate [fj|, the radiation from the 
plate, given in a plane wave basis, must be transformed 
to the spherical basis 18(, reflected by the sphere, trans- 
formed back and so on. For simplicity, we focus on a plate 
held at a finite temperature T p ^ 0, while the sphere and 
environment are at zero temperature ( T s = T env = 0). 
This suffices to describe also situations with T s ^ as 
the transfer vanishes for T p = T s , and hence for T s ^ 0, 
we subtract our result evaluated at T s . We express the 
correlations in Eq. ((3]) in a plane waves basis, and H s , the 
energy absorbed by the sphere, is obtained by integrating 
S in Eq. ^ over two infinite parallel planes enclosing the 
sphere and separating it from the plate. 

Figure [2] shows the results for the heat transfer from 
a Si02 plate at room temperature to a SiC>2 sphere of 
R = 5/iin at zero temperature, with surface-to-surface 
separation d, normalized by the Stefan-Boltzmann law 
H s = crT' l 2irR 2 (only half of the sphere is exposed to 
the plate). For large d, H s is roughly 0.5 in these units, 
whereas for d — >• 0, H s diverges due to the increased 
tunneling of evanescent waves, eventually exceeding the 
Stefan-Boltzmann value. The figure shows the numerical 
solution of Eq. ^ together with a one reflection approx- 
imation, where we set O ai( g = (1 — GoTg) in Eq. ([5]), 
neglecting higher order reflections between sphere and 
plate. We see that the two curves approach each other 
for large d, as most rays are scattered outward and will 
not hit the sphere a second time. The reflection ex- 
pansion is hence helpful for getting analytical results for 
d R. Our numerical solution involves an expansion 
in spherical multipoles: For R/Xt large or d/R small, 
increasingly more multipoles arc needed. In practice, we 
restrict to a maximal multipolc order of l max = 20, for 
accurate results up to d > R/2. Since closer separations 
are also interesting and relevant experimentally, but diffi- 
cult numerically, we demonstrate in the inset of Fig.[2]thc 
approach to a proximity transfer approximation (PTA), 
equivalent to the proximity force approximation (PEA) 
used in Casimir physics, 



lim HJd) 

d/R^O 



2ttR 



d+R 



S pp {s)ds, 



(8) 



other descriptions have been offered [22 1 



where S pp (d) is the Poynting vector for parallel plates at 
separation d. We identify the divergent terms as d — s- 
for both the sphere-plate and plate-plate configurations 
(the E modes originating from evanescent waves), and 
evaluate their ratio in the one reflection approximation 
(allowing us to use l max = 200). As demonstrated in 
Fig. [21 this ratio approaches unity for d — > 0, suggesting 
T of the sphere approaches T of the plate in the PTA- 
sense. From this, we anticipate that multiple applications 
of these matrices (leading to the full solution) will also 
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FIG. 2. Heat transfer rate (in units of Stefan-Boltzmann's 
law) from a room temperature plate to a sphere at T = 
of radius R — 5/j.m (both Si02), as a function of separation. 
The horizontal bar shows the limit of d — >• oo. The left inset 
shows the approach to PTA for the divergent terms in a one 
reflection approximation. The right inset shows the result at 
large separation as function of R. 



approach the ratio unity (independent of the accuracy of 
the one reflection approximation as d — > 0). We investi- 
gated different R and confirmed that PTA in Eq. fjSJ is 
valid in general, with H s (d) oc d _1 as d — > 0. While a 
similar point is discussed in Ref. @ for the case of two 
spheres and used in experimental studies @, 0] , to our 
knowledge the validity of PTA was not quantitatively 
demonstrated previously. It is not obvious as it implies 
that the ratios R/Xt and R/5 are irrelevant as d — >• 0. 

For d — >• oo, H s approaches a constant, which can be 
obtained by considering the d independent part of the 
plate radiation, and using the one reflection term. The 
result, as shown in the right inset of Fig. [21 is quite similar 
to the behavior of a single sphere in Fig. Q] For small R, 
H s is proportional to the volume of the sphere, for similar 
reasons as discussed before. In this limit, H s is given by 
(with magnetic permeability of the sphere [i s and Fresnel 
reflection coefficients r E and r M of the plate for angle 9) , 



lira iJ, = 



cR J 



d>A T >-R 167T 3 J 

r-Tr/2 



dujaxioj) Im 



Ms - 1 



1 



dO sin 9 



P=E.M 



fi s + 2 e s + 
(l-\r p (9,u>)\ 2 ). (9) 



For R ^> At, we may expect the result to approach a 
classical limit, given by crT 4 e 2 27ri? 2 , with e from Fig. 1. 
While the data points come close to this value, one does 
not expect exact approach Q , in contrast to Fig. [1] be- 
cause the Fresnel coefficients depend on the angle of in- 
cidence. If we additionally let (e p ,e s ) — >• 1, H s will ap- 
proach the classical limit since the Stefan-Boltzmann law 



applies to all convex black bodies. 

While we highlighted applications to simple shapes, 
the formalism presented here is more general, and com- 
bined with a numerical scheme for the computation of 
scattering matrices 23j can deal with collections of ob- 
jects at different temperatures. Indeed, such a formal- 
ism is needed to properly deal with near field effects in 
device and fabrication at the micron scale. The formal- 
ism yields also Casimir forces between objects at differ- 
ent temperatures- examples of which we leave for future 
work. We note, however, that in the final stages of this 
project we became aware of two independe nt, p artly re- 
lated, studies of non-equilibrium effects 24 125|. 
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